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We review the effects of electron-electron interactions on the ground-state spin and the transport 
properties of ultra-small chaotic metallic grains. Our studies are based on an effective Hamiltonian 
that combines a superconducting BCS-like term and a ferromagnetic Stoner-like term. Such terms 
originate in pairing and spin exchange correlations, respectively. This description is valid in the 
limit of a large dimensionless Thouless conductance. We present the ground-state phase diagram in 
the fluctuation-dominated regime where the single-particle mean level spacing is comparable to the 
bulk BCS pairing gap. This phase diagram contains a regime in which pairing and spin exchange 
correlations coexist in the ground-state wave function. We discuss the calculation of the tunneling 
conductance for an almost-isolated grain in the Coulomb-blockade regime, and present measurable 
signatures of the competition between superconductivity and ferromagnetism in the mesoscopic 
fluctuations of the conductance. 



INTRODUCTION 

Superconductivity and ferromagnetism compete with 
each other. Pairing correlations lead to Cooper pairs of 
electrons with opposite spins and thus tend to minimize 
the total spin of the grain, while ferromagnetic correla- 
tions tend to maximize the total spin. 

Nevertheless, it is well known that superconducting 
and ferromagnetic order can be present simultaneously 
in bulk systems when ferromagnetism is caused by local- 
ized paramagnetic impurities 

0, HHH S- Recently, 
it was observed that both states of matter can coexist 
in high-Tc superconductors 0, 0] and in heavy fermion 
systems [H 0, [l(| even when the electrons that are re- 
sponsible for superconductivity and ferromagnetism are 
the same. This surprising observation led to the search 
of new theoretical models that can describe this coexis- 
tence. 

In ultra-small metallic grains, in which the bulk pair- 
ing gap A is comparable to the single-particle mean level 
spacing 5, a coexistence regime of superconductivity and 
ferromagnetism was predicted [H El El- The ground 
state of the grain is described by a state where a few 
single-particle levels around the Fermi energy are singly 
occupied while all other electron are paired. This coexis- 
tence regime is characterized by spin jumps and its size 
can be tuned by an external Zeeman held. 

However, it is difficult to measure the ground-state spin 
of a grain, and a more directly measurable qua ntity is the 
tunneling conductance through the grain UJ. In addi- 
tion, one has to take into account the mesoscopic fluc- 
tuations that are typical for chaotic grains [lij]. Effects 
of exchange correlations on the conductance statistics in 
quantum dots, in which pairing correlations are absent, 
were studied in Ref. 

Eft 

In Ref. we identified signa- 



tures of the coexistence of pairing and exchange correla- 
tions in the mesoscopic fluctuations of the conductance 
through a metallic grain that is weakly coupled to leads. 

The fabrication and control of nano-sizc metallic de- 
vices is a challenging task. The first conductance mea- 
surements in ultra-small metallic grains were carried out 
in the mid-nineties [3, [l9|, [2(|. The grains were pro- 
duced by breaking nanowires and their size was diffi- 
cult to control. Coulomb blockade, discrete levels and 
pairing effects were observed in a single grain by mea- 
suring the tunneling conductance 14[. During the last 
decade numerous technological advances led to an in- 
crease in control and tunability of ultra-small metallic 
grains. Break junction techniques 21] and electromigra- 
tion were used for gating and establishing precise 
contact between leads and grain. A particularly impor- 
tant recent development has been the use of monolayers 
of organic molecules as tunnel barriers, enabling control 
of the size and shape of the grain [231 ]. New materi- 
als have been tested as well. Cobalt nanoparticles were 
used to investigate the effect of ferromagnetism 24, 25|. 
Spin-orbit coupling and non-equilibrium excitations were 
studied in gold grains [22I liij . [2r3 |. The recent discov- 
ery of superconductivity in doped silicon at atmospheric 
pressure and critical temperatures of a few hundred mil- 
likelvin [23] might further facilitate the development of 
mesoscopic superconducting devices. However, the com- 
petition between superconductivity and ferromagnetism 
has not been investigated so far. 

Here we review the effects of electron-electron interac- 
tions on the ground-state spin and transport properties of 
ultra-small metallic grains. Our analysis is based on an 
effective Hamiltonian for chaotic or disordered systems 
that combines a BCS-like pairing term and a Stoner-like 
spin exchange term. This so-called universal Hamilto- 
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nian [28|, [29[ is valid in the limit of a large Thouless con- 
ductance. This universal Hamiltonian and its solution 
are described in Sec. . In Sec. we present the phase 
diagram of the ground-state spin and discuss a regime in 
which superconductivity and ferromagnetism coexist. In 
Sec. we review the mesoscopic fluctuations of the tun- 
neling conductance through an almost-isolated metallic 
grain. In particular, we discuss signatures of the coex- 
istence of pairing and exchange correlations in the con- 
ductance peak height and peak spacing statistics. We 
also propose specific materials for which such mesoscopic 
coexistence might be observed experimentally. 

MODEL 

We consider a chaotic metallic grain with a large di- 
mensionless Thouless conductance. The low-energy ex- 
citations of such a grain are described by an effective 
universal Hamiltonian [28|, [2i| 

H = e * c L c ^ - Gp ^P - , (!) 

ka 

where c kcr is the creation operator for an electron in 
the spin-degenerate single-particle level efc with spin up 
(cr = +) or spin down (a — — ). The first term on the 
r.h.s. of {T]) describes the single-particle Hamiltonian of 
an electron in the grain (i.e., kinetic energy plus confining 
single-particle potential). The second term is the pairing 
interaction with coupling constant G and the pair cre- 
ation operator pt — Si c !+ C l-- The third term is an 

exchange interaction where S = Ylkaa' c'ka T oo' c ^ 1 x& * ne 
total spin operator [ji are Pauli matrices) and J s is the 
exchange coupling constant. Estimated values of J s for a 
variety of materials were tabulated in Ref. [30( ■ In Eq. (fl]) 
we have omitted the charging energy term e 2 N 2 /2C (C 
is the capacitance of the grain and N is the number of 
electrons) since it is constant for a grain with a fixed 
number of electrons. 

The universal Hamiltonian ([1]) conserves the total spin 
of the grain, i.e., [H,S] = 0. Consequently, each eigen- 
state has a well-defined total spin S and spin projec- 
tion M. The pairing interaction scatters pairs of spin 
up/down electrons from doubly occupied to empty or- 
bitals. Therefore the pairing interaction does not affect 
the singly occupied levels. This is known as the blocking 
effect and the singly occupied levels are also referred to as 
'blocked levels'. On the other hand, these singly occupied 
levels are the only levels that contribute to the exchange 
interaction. Thus, each eigenstate of (JTJ) factorizes into 
two parts. The first part \IA) is a zero spin eigenstate 
of the reduced BCS Hamiltonian J2ka e fc c iL c fc<x — GP^P, 
and is described as a superposition of Slater determinants 
that are constructed from the subset U of empty and dou- 
bly occupied levels. The second part of the eigenstate, 



\B, 7, S, M), is obtained by coupling the set of singly oc- 
cupied levels B, each carrying spin 1/2, to total spin S 
and spin projection M [161 l3l| . Here, 7 denotes a set 
of quantum numbers distinguishing between eigenstates 
with the same spin and singly occupied levels. For a given 
set B of b singly occupied levels, the allowed spin values 
vary between S — (S — 1/2) for even (odd) number of 
electrons and S = 6/2. Each of these spin values has a 
degeneracy of 

The complete set of eigenstates is then given by 

\i) = \U,B n ,S,M). (3) 

The reduced pairing Hamiltonian is characterized by 
a coupling constant G. However, the physical parameter 
that determines the low-energy spectrum of the grain (for 
J s = 0) is A/ 5, where A is the bulk pairing gap and 5 
the single-particle mean level spacing. We can truncate 
the total number of levels from N to N r < N , and 
renormalize G such that the low-energy spectrum of the 
grain remains approximately the same. For a picketfence 
spectrum, the renormalized coupling constant is given by 



arcsmh [-^jf- ) 

The exchange interaction only affects the singly occu- 
pied levels, and we expect the renormalization @ to hold 
as long as the number of singly occupied levels is small 
compared with N r . In practice, we study the relevant 
observable as a function of truncated bandwidth N r and 
make sure that its value has converged for the largest 
bandwidth N r we can calculate. 

GROUND-STATE PHASE DIAGRAM 

In this section we study the ground-state spin of the 
grain as a function of J s /S and A/S. We find the lowest 
energy E(S) for a given spin S and then minimize E(S) 
with respect to S. The ground-state spin of the grain is 
determined by the competition between various terms in 
the Hamiltonian (JTJ) - The one-body part and pairing in- 
teraction favor minimal spin, while exchange interaction 
favors a maximally polarized state. 

The ground-state phase diagram in the A/S-J s /5 
plane of a grain with a picketfence single-particle spec- 
trum is shown in Fig.[T] We find three different phases: a 
superconducting phase (SC) where the number of pairs is 
maximal and S is minimal, a ferromagnetic phase (FM) 
where the system is fully polarized S = N/2, and an 
intermediate regime (SC-FM) where exchange and pair- 
ing correlations coexist. The ground-state wave function 
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FIG. 1: Ground-state phase diagram of a grain with picket- 
fence single-particle spectrum in the J s /5-A/S plane for an 
even (left panel) and odd (right panel) number of electrons. 
Numbers denote the spin values in the corresponding sectors. 
We observe an intermediate regime (SC-FM) in which the 
ground state is partly polarized and partly paired. Taken 
from Ref. [131 ]. 



with spin S in the coexistence regime is described by 
b = 2S singly occupied levels closest to the Fermi energy 
while the remaining electrons are paired. 

The coexistence regime is bounded by two critical val- 
ues Js /5 and Jg/5 of the exchange interaction that 
are function of A/d. The lower value Js/S is a mono- 
tonically increasing function of A /S (stronger exchange is 

required to polarize a grain with stronger pairing correla- 
te^ 

tions), while the higher value Js /6 depends only weakly 
on A/ 5. 

It is interesting to follow the dependence of the ground- 
state spin as a function of the exchange coupling constant 
J s /5 for a fixed value of A/S. In the absence of pairing 
(A = 0), this dependence follows a stepwise behavior 
known as the mesoscopic Stoner staircase (28|, where a 
transition from spin S to spin S+l occurs for an exchange 
coupling of 



S 



2S+1 
2S + 2 



at A = . 



(5) 



The first step occurs at J s /S = 0.5 (where the ground- 
state spin increases from S = to S — 1) and continues 
up to J s — 0.75 (where the S = 1 to S = 2 transition 
takes place). In the presence of pairing, the first step is 
shifted to higher values of the exchange and the Stoner 
staircase is compressed. For A/ 5 < 0.6, all steps have a 
height of AS = 1. However, for 0.6 < A/S < 0.8, the 
first step has a height of AS* = 2, describing a spin jump 
from S = to S = 2. This first step starts at J s /S « 0.87 
and ends at J s /8 « 0.9. The height of the first-step spin 
jump increases at larger values of A/ 5. All subsequent 
steps arc of height one. 

Spin jumps also occur when superconductivity in 
metallic grains breaks down in the presence of a suffi- 



ciently large external Zeeman field 32J. This "softened" 
first-order phase transition from a superconductor to a 
paramagnet was explained qualitatively using a finite- 
spin BCS approximation. 

In the presence of exchange correlations, spin jumps 
are predicted to occur at J s /S > 0.87. Such exchange 
coupling values are significantly larger than the exchange 
coupling values of most metals (see Fig. 9 in Ref. [30l|). 
Moreover, the exchange strength is an intrinsic mate- 
rial property and is difficult to tune experimentally. In 
Ref. [13| we have shown that the coexistence regime can 
be tuned to experimentally accessible values of the ex- 
change interaction by applying an external Zeeman field. 



CONDUCTANCE 

The determination of the ground-state spin of a grain 
is a difficult measurement. It is then desirable to identify 
signatures of coexistence of superconductivity and ferro- 
magnetism in a quantity that is directly measurable, e.g., 
the conductance. Furthermore, the universal Hamilto- 
nian JT]) is only valid for chaotic (or disordered) grains, 
in which mesoscopic fluctuations are generic. Therefore, 
in order to make quantitative predictions for the exper- 
iment it is necessary to include the effect of mesocopic 
fluctuations. In this section we discuss the mesoscopic 
fluctuations of the tunneling conductance for an almost- 
isolated metallic grain in the Coulomb blockade regime. 
We find signatures of the coexistence of pairing and ex- 
change correlations in the conductance statistics. Since 
the tunneling conductance can be measured in a single- 
electron transistor that uses the metallic grain as an is- 
land, our results are directly relevant for the experiment. 

We consider grains that are weakly coupled to external 
leads. In the regime of sequential tunneling S, T s> T (T 
is a typical tunneling width). Assuming the charging en- 
ergy to be much larger than temperature (e 2 /2C s> T), 
the conductance displays a series of sharp peaks as a 
function of gate voltage. The A-th conductance peak de- 
scribes a tunneling event in which the number of electrons 
in the dot changes from A to N + 1. and is determined 
by the many-body energies and transition rates between 
eigenstates of the A and A + 1 electrons. The conduc- 
tance peak height and peak position are determined by 
solving a system of rate equations [33j]. 

Here we present results for the conductance peak spac- 
ing and peak height statistics for an experimentally acces- 
sible temperature of T = 0.1 S [20(. In the absence of an 
external magnetic field, the single-particle Hamiltonian is 
described by the Gaussian orthogonal ensemble (GOE) 
of random matrices. For each random matrix realiza- 
tion of the one-body Hamiltonian, we use the Lanczos 
method to find the five lowest eigenstates of the many- 
body Hamiltonian ([T]). The calculations are carried out 
for a truncated bandwidth A^. r = 8 and electron numbers 




FIG. 2: Peak spacing distributions at T — 0.1 S for several 
values of A/S. Results are shown for both J s = (dashed, 
grey-filled histograms) and J a — 0.6 8 (solid histograms). For 
A = we also compare with the analytic result [341 ] at T <C 5 
and J s = (solid line) . The bimodality of the distribution at 
A = J s = is suppressed at finite exchange (J s = 0.6 S) but 
is restored for A/S = 0.5. Taken from Ref. [171 ]. 
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FIG. 3: Peak height distributions at T — 0.1 S for several 
values of A/S. We show results for J s = (dashed, grey- 
filled histograms) and J s = 0.6 5 (solid histograms). The 
solid line is the analytic distribution at T -C 5 and J s = [34| . 
Since the conductance peak height (? max fluctuates over several 
order of magnitude, we show the distributions as a function 
of ln((/ max /<7 max ) where <? max is the average conductance peak 
height. Taken from Ref. \vft. 



N = 16, 17, 18 and 19. Using the many-body energies 
and wave functions, we calculate the tunneling matrix 
elements for the corresponding three tunneling events. 
We then solve the system of rate equations and deter- 
mine the conductance as a function of gate voltage. The 
peak position and height are determined numerically. To 
ensure good statistics, the above procedure is repeated 
for 4000 realizations of the one-body Hamiltonian. 

Peak spacing statistics 

The peak spacing distribution is shown in Fig. [51 where 
the spacing is measured relative the constant charging 
energy. For A = J s = and at low temperatures, this 
distribution is bimodal because of the spin degeneracy of 
the single-particle levels [l5j]. The exchange interaction 
induces mesoscopic spin fluctuations and suppresses this 
bimodality (see top left panel of Fig. [5]) . This is known 
from the case of semiconductor quantum dots. 

Pairing correlations can restore bimodality. For a mod- 
erate exchange value of J s — 0.3 5, bimodality is restored 
for a relatively weak pairing strength of A = 0.25 5. For 
J s = 0.6 S, this bimodality is suppressed but reappears at 
A/S = 0.5 (see bottom left panel of Fig.HJ). The left part 
of the peak spacing distribution describes even-odd-even 
(E-O-E) tunneling events (the parity refers to the number 
of electrons), and its right part describes odd-even-odd 
(O-E-O) transitions. 



The bimodality of the peak spacing distribution in 
the presence of strong pairing correlations can be un- 
derstood qualitatively in the T — fixed-spin BCS ap- 
proximation [l3j]. For an E-O-E transition, the first con- 
ductance peak corresponds to the blocking of an addi- 
tional single-particle level, while the second conductance 
peak corresponds to the removal of this blocked level 
by creating an additional Cooper pair. This leads to 
the estimate A EOE w -2A + |j. In a O-E-0 tunnel- 
ing sequence, these two events are reversed and we find 
A^ EO « 2 A - | J. The contribution of the exchange 
interaction in these estimates is straightforward because, 
in the limit of strong pairing, the ground-state spin is 
always S = (S = 1/2) for an even (odd) number elec- 
trons. The difference of these two peak spacing values 
is 



SA 2 = A° EO - A E w 4A - 3 J , A > S , (6) 



and bimodality becomes more pronounced when A/ 8 in- 
creases. Since the exchange interaction strength for most 
metals is smaller or comparable to J s ~ 0.6 6, exchange 
correlations are insufficient to suppress the bimodality in 
the presence of strong pairing correlations. 
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Peak height statistics 

The peak height distribution is shown in Fig. [3J For 
A = J s = 0, this distribution is known analytically at 
T <C 5 and, in the absence of an external magnetic 

field (GOE statistics), diverges at small values of the con- 
ductance. Finite temperature and exchange interaction 
have a similar effect on the peak height distribution; they 
both reduce the occurrence of small conductance values. 
While thermal fluctuations open an energy window in 
which states become available for tunneling and thus 
can contribute to the conductance, the exchange inter- 
action increases the many-body density of states around 
the Fermi energy and makes otherwise high-lying non- 
zero spin states available for tunneling. This effect is 
shown in Fig. [3] For A = we clearly observe that 
a finite exchange interaction suppresses the peak height 
distribution at small conductance values. The same be- 
havior was observed for the Gaussian unitary ensemble 
(GUE) ensemble in semiconductor quantum dots, where 
the pairing interaction can be ignored and a closed so- 
lution for the conductance is available [3]. There, the 
suppression of probability at small conductance values 
by the inclusion of exchange interaction leads to better 
agreement with the experimental results. Another signa- 
ture of exchange correlations is the suppression of peak 
height fluctuations as described by o{g- m&yL )/g max (the ra- 
tio between is the standard deviation and average value 
of the peak height) [ij}- 

However, it is important to note that at very low tem- 
peratures, e.g., T = 0.01 5, the small conductance values 
are no longer suppressed by an exchange interaction. For 
such very low temperatures, the increase in the density of 
states near the Fermi energy by the exchange interaction 
does not affect the conductance since only the ground 
state contributes significantly to transport. 

The pairing interaction leads to an excitation gap that 
pushes states with large spin to higher energies. Thus, 
already for A = 0.5 6 we observe less suppression of the 
small conductance values by an exchange of J = 0.6 6. 
When the pairing interaction is strong enough to suppress 
all spin polarization of low-lying states (e.g., A > 1.0 6 
for J s < 0.6(5), the peak height distribution becomes 
essentially independent of the exchange interaction and 
diverges again for small conductance values. 

For A/6 = 0.5 and J s /6 — 0.6, we observe a signa- 
ture of pairing correlations in the peak spacing distri- 
bution (bimodality) and a signature of exchange corre- 
lations in the peak height distribution (suppression of 
the probability of small conductance values). We can 
interpret these results to describe the mesoscopic coex- 
istence of pairing and exchange correlations. A candi- 
date for this mesoscopic coexistence is vanadium, which 
has J s /6 ~ 0.57 — 0.63 [3(| and is superconducting in 
the bulk. Another candidate is platinum, which has 



J s /6 ~ 0.59 — 0.72 and is superconducting in granular 
form. 



CONCLUSION 

We reviewed the competition between superconductiv- 
ity and ferromagnetism in ultra-small metallic grain. In 
particular, we presented the ground-state phase diagram 
in the J s /6 — A/5 plane, and discussed a coexistence 
regime of superconductivity and ferromagnetism. This 
regime is characterized by spin jumps that are greater 
than unity. We also discussed the transport properties 
of the grain in its Coulomb blockade regime of weak cou- 
pling to leads and described the statistics of the con- 
ductance peaks in the presence of both pairing and ex- 
change correlations. Of particular interest is a regime 
in which pairing and ferromagnetic correlations coexist. 
Such a regime is defined by the simultaneous occurrence 
of bimodality in the peak spacing distribution (caused 
by pairing correlations) and the suppression of the peak 
height distribution at small conductance values (caused 
by ferromagnetic correlations). 

The exchange interaction strength is a material con- 
stant and might be difficult to tune experimentally. Al- 
ternatively, the coexistence regime can be controlled by 
an external Zeeman field while measurements are carried 
out at a fixed exchange interaction strength. 
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grant No. DE-FG-0291-ER-40608. Computational cycles 
were provided by the NERSC high performance comput- 
ing facility at LBL and the Yale Bulldog cluster. 
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